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Abstract 



Simple, self-similar, analytic solutions of 1 + 1 dimensional relativistic hydrody- 
namics are presented, generalizing the Hwa - Bjorken boost-invariant solution to 
inhomogeneous rapidity distributions. These solutions are generalized also to 1 + 
3 dimensional, cylindrically symmetric firetubes, corresponding to central collisions 
of heavy ions at relativistic bombarding energies. 
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1 Introduction 

Analytic solution of the equations of relativistic hydrodynamics is a difficult 
task because the equations are non-linear partial differential equations, that 
are rather complicated to handle not only analytically but also numerically. 
However, relativistic hydrodynamics has various applications, including the 
calculations of single-particle spectra and two-particle correlations in relativis- 
tic heavy ion collisions, see ref [1] . More recently, there has been an increasing 
interest in applications of relativistic hydrodynamics in Au+Au collisions at 
RHIC both at y/s = 130 AGeV and \fs = 200AGeV bombarding energies, 
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predictions were made for the coming LHC experiments [2,3,4]. The hydro- 
dynamical analysis can also be extended to the study of these processes on 
event-by-event basis [5,6]. However, most works in hydrodynamics are numer- 
ical so not always transparent. 

In this sense, exact solutions would be useful, but are rarely found due to 
the highly non-linear nature of relativistic hydrodynamics. Khalatnikov's one- 
dimensional analytical solution [7] to Landau's hydrodynamic model [8] gave 
rise to a new approach in high energy physics. The boost-invariant solu- 
tion [9] was found later by R. C. Hwa and other authors. It has been fre- 
quently utilized as the basis for estimations of initial energy densities in ultra- 
relativistic nucleus-nucleus collisions [10]. Due to this famous application this 
boost-invariant solution is frequently called as Bjorken's solution, although as 
far as we know it was first described by R. Hwa in ref. [9]. Perhaps it should 
be called the Hwa-Bjorken solution, which name we shall use hereafter. 

Recently, Biro has found self-similar exact solutions of relativistic hydrody- 
namics for cylindrically expanding systems [11,12]. However, his solutions are 
valid only when the pressure is independent of space and time, as e.g. in the 
case of a rehadronization phase transition in the middle of a relativistic heavy 
ion collision. 

Here we present an analytic approach, which goes back to the data-motivated 
exact analytic solution of non-relativistic hydrodynamics found by Zimanyi, 
Bondorf and Garpman (ZBG) in 1978 for low energy heavy ion collisions 
with spherical symmetry [13]. This solution has been extended to the case of 
elliptic symmetry by Zimanyi and collaborators in ref. [14]. In [15,16] a Gaus- 
sian parameterization has been introduced to describe the mass dependence 
of the effective temperature and the radius parameters of the two-particle 
Bose-Einstein correlation functions in high energy heavy ion collisions. Later 
it has been realized that this phenomenological parameterization of data cor- 
responds to an exact, Gaussian solution of non-relativistic hydrodynamics 
with spherical symmetry [17]. The spherically symmetric self-similar solutions 
of non-relativistic hydrodynamics were obtained in a general manner in [18], 
that included an arbitrary scaling function for the temperature profile, and 
expressed the density distribution in terms of the temperature profile func- 
tion. The ZBG solution and the Gaussian solution of [17] are recovered from 
the general solution of [18] as special cases, corresponding to different scaling 
functions of the temperature profile. The Gaussian solution has been gener- 
alized to ellipsoidal expansions in [19], that provides analytic insight into the 
physics of non-central heavy ion collisions [20]. 

Our approach corresponds to a generalization of these recently obtained an- 
alytic solutions [17,18,20,21] of non-relativistic fireball hydrodynamics to the 
case of relativistic longitudinal and transverse flows. In particular, an analytic 
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approach, the Buda-Lund (BL) model has been developed to parameterize 
the single particle spectra and the two-particle Bose-Einstein correlations in 
high-energy heavy-ion physics in terms of hydrodynamically expanding, cylin- 
drically symmetric sources [22]. Here we attempt to find a family of exact 
solutions of relativistic hydrodynamics that may include the BL model as a 
particular limiting case. It turns out that in the simplest case our result cor- 
responds to the Cracow hydrodynamic parametization, which is successfull in 
describing single particle spectra of Au+Au collisions at = 130 and 200 
AGeV at RHIC[23,24,25]. 



2 The equations of relativistic hydrodynamics 



We solve the relativistic continuity and energy-momentum conservation equa- 
tion: 

d>o = o, (i) 

d v T» v = . (2) 

Here n = n(t, r) is the number density, the four-velocity is denoted by w M = 
■u M (t, r) = 7(1, v), normalized to u^u^ = 7 2 (1 — v 2 ) = 1, and the energy- 
momentum tensor is denoted by T^ u . We assume perfect fluid, 

= (e + p)uV - pg^, (3) 

where e stands for the relativistic energy density and p denotes the pressure. 

We close this set of relativistic hydrodynamical equations with the equations 
of state. We assume a gas containing massive conserved quanta, 

e = mil + Kp , (4) 
p = nT. (5) 

The equations of state have two free parameters, m and k. Non-relativistic 
hydrodynamics of ideal gases corresponds to the limiting case of m > T, 
v 2 < 1 and k — 3/2 . Relativistic hydrodynamics for massless particles and a 
constant speed of sound c 2 corresponds to the case of m = and c 2 = 1/k. 

The energy-momentum conservation equations can be projected into a com- 
ponent parallel to and components orthogonal to u^, which are respectively 
the relativistic energy and Euler equations: 

u^e+(e + P )3/ = 0, (6) 
u^d^p + (e + p)u> 1 d^u u - d v p = . (7) 
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Based on general thermodynamical considerations, one can show that the 
expansion is adiabatic: 

dX) = 0, (8) 

where a is the entrophy density. This relation holds for perfect fluids, inde- 
pendently of the equations of state. 

With the help of the equations of state and the continuity equation, the energy 
equation can be rewritten as an equation for the temperature, 

u^T + -Td^ = . (9) 



We solve 5 independent equations, the continuity, the (3 spatial components 
of) relativistic Euler, and the temperature equation, eqs. (1,7,9). The equa- 
tions of state, eq. (4,5) close this system of equations in terms of 5 variables, 
n, T and v = (v x ,v y ,v z ). 



3 Self-similarity 



We look for solutions which generalize the usual similarity flow, in which the 
flow pattern is unchanged with time if the scales of length X(t), Y(t), Z(t) 
along three orthogonal directions vary appropriately, namely, we consider 

\X(t) x 'Y(t) v 'Z(t) z )' 1 J 

where = (t, r x ,r y ,r z ) and the dot indicates the time derivative. As for the 
thermodynamic quantities such as n(x M ), T(x M ),p(a; M ), . . . , we search solutions 
of the form 

f(^) = fo(^) a F(s), (11) 



V. 

where the volume parameter V = XYZ, a is an appropriate exponent and 
F(s) is an arbitrary fuction of the scaling variable defined by 

s= # + # + #- (12) 



These are Hubble type of flows, but the thermodynamic quantities may contain 
arbitrary functions depending on the the scale parameter s and also, at least 
in principle, the scale parameters X(t),Y(t) and Z(t) may be different in 
the principal directions. Their derivatives, X(t), Y(t) and Z(t) correspond to 
(direction and time dependent, generalized) Hubble constants. 
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In heavy-ion collisions, the well known boost-invariant solution [9] is often 
utilized to discuss several properties of data. However, this solution has some 
shortcomings: i) it is scale invariant, having a flat rapidity distribution, cor- 
responding to the extreme relativistic collisions; ii) it contains no transverse 
flow. In the present paper, we apply the strategy described above first to 
1 + 1 dimensional (time + longitudinal coordinate) case and obtain a class 
of solutions which are able to describe inhomogeneous rapidity distributions, 
overcoming the first shortcoming mentioned above. Then, in section 5, we con- 
sider the case of cylindrically symmetric case, trying to overcome the second 
shortcoming. 



4 Simple 1+1 dimensional solutions 



In this section, we solve the 1+1 dimensional problem. Hence = (t, r z ), = 
(E, k z ) throughout this section. The metric tensor is g^ v = = diag(l, — 1) 
and rr M = (t,—r z ). We solve 3 independent equations, the continuity, the 
temperature equation and the z component of the Euler equation (1,7,9). The 
equations, (4,5) close this system of equations in terms of 3 variables, n, T 
and v z . 

We look for flows that scale in the z direction. The scaling variable, eq.(12), 
in this case is defined as 

(13) 



and the longitudinal velocity 

u *(*> r *)=fH r *> (i4) 

where Z = dZ(t)/dt . In the relativistic notation, this form is equivalent to 

= (coshC, sinhC), (15) 
tanhC = ^7^r z , or coshC = —. 1 =7- (16) 

Z W yjl - Z 2 S 

Note that from eq. (16) it is obvious that this solution can be defined only in a 
bounded longitudinal coordinate region, because at any time \r z \ < Z(t)/Z(t) 
has to be satisfied. Using this ansatz, we find that the continuity equation is 
solved by the form 

n(t,r 2 )=n §J^( S ), (17) 

where Q(s) is an arbitrary non- negative function of the scaling variable s and 
n and Z are normalization constants. We use the convention Z = Z(t ) and 
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no = n(to, 0) which implies that Q(s = 0) = 1. The temperature equation, (9) 
is solved by the following form: 




T(t,r z )=T F(s). (18) 



The constants of normalization are chosen such that T = T(t , 0) and .F(O) = 
I. Here again, we find that the solution is independent of the form of the 
function J-'(s). From the positivity of the temperature distribution it follows 
that F(s) > 0. 

Using the ansatz for the flow profile and the solution for the density and the 
temperature, the relativistic Euler equation reduces to a complicated non- 
linear equation that contains Z, Z and Z and s. Taking this equation at s = 
we express Z as a function of Z and Z. Substituting this back to the Euler 
equation we obtain an equation for Z, Z and s. In particular, for the m = 
case, Z cancels out and this reduces to a second order polynomial equation 
for Z 2 , which has only one positive root. The form of the solution in this case 
(m = 0) is Z 2 (t) = F(s). Observing that the function F depends only on the 
scaling variable s, while Z depends only on the time variable t, we conclude 
that the only solution of this equation should be a constant Z = Zq. Now 
we choose the origin of the time axis such that Z(t — 0) = without loss of 
generality. The solutions can be cast in a relatively simple form by introducing 
the longitudinal proper time r and the space-time rapidity rj, 



V = hog( t -±^). (20) 



t 2 -r 2 , (19) 

This implies that Z(t) = Z t, v z — y = tanh^ and ( = rj. Thus the solution 
for the flow velocity field corresponds to the flow field of the boost-invariant 
solution. However, in the boost-invariant solution the temperature distribu- 
tion was independent of the 77 variable, while in our case the density and the 
temperature distributions can be both i] dependent, or in other words, our 
solutions are scale dependent. The scale is defined by the parameter Z , in the 
longitudinal direction. 

This special form of the solution for the flow velocity field implies that Z = 
0. This equation implies that there is no pressure gradient and there is no 
acceleration in this class of self-similar solutions, similarly to the case of boost- 
invariant solution. The Euler equation is reduced to the following requirement: 



(d z + r jd t ) 



/+ \ (1+1 Ik) 

(£) {i-ztst+^g{s)Hs) 



(21) 



This equation is solved by the trivial Q{s)F{s) = as well as by the non-trivial 
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solution of 

g(s)^(s) = (l-Z 2 s)-^\ (22) 

which is indeed only a function of s as Z is a constant of time. With this 
form, the Euler equation is satisfied. This solution implies that the scaling 
profile functions for the temperature and the density distribution are not in- 
dependent. As the constraint is given only for their product, one of them can 
be still chosen in an arbitrary manner. 

It is worthwhile to introduce new forms of the scaling functions. Let us define 

T{s) = Hs){l-Zlsf'^ (23) 
V(s) = G(s)(l - Z 2 s) . (24) 

Then the constraint Eq. (22) can be cast to the simplest form of 

V(s)T(s) = 1. (25) 

Let us summarize our new family of solutions of the 1+1 dimensional rela- 
tivistic hydrodynamics by substituting the results in the density, temperature 
and pressure profiles. We obtain 

T 

v z — — — tanh?7, (26) 



r\ tanh 2 r\ 



Zot 2 Zq 
to- 



(27) 



n = n -V(s), (28) 

T 

,f \ 1+1/k 

P = Po ) , (29) 

/f \ 1/k i 

T^n(±) ^, ,30) 

where po — n-oTo- Thus we have generated a new family of exact solutions 
of relativistic hydrodynamics: a new hydrodynamical solution is assigned to 
each non-negative function V(s). It can be checked that the above solutions 
are valid also for massive particles, the form of the solution is independent of 
the value of the mass m. The form of solutions depends parametrically on k, 
that characterizes the equation of state. 



4-1 Analysis of the solutions 



The pressure and the flow profiles of the above 1+1 dimensional relativistic 
hydro solution are the same as in the boost-invariant solution. In the case of 
V(s) = 1, we recover the Hwa-Bjorken boost-invariant solution of refs. [9,10]. 
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In this limiting case, the pressure, the density and the temperature profiles 
depend only on the longitudinal proper time r. 

In the general case, our solution contains a characteristic scale defining pa- 
rameter in the longitudinal direction, Z , and an arbitrary scaling function 
V(s). Thus we have an infinitely rich new family of solutions. Let us try to 
determine the physical meaning of the scaling function V(s). 

In order to do this we evaluate the single particle spectra corresponding to 
the new solutions. Here we neglect any possible dynamics in the transverse 
directions, as usual in case of applications of the boost-invariant solution. 
The four- velocity field of our solutions thus becomes w M = (cosh 77, 0, 0, sinhr/). 
The four-momentum of the observed particles with mass m is denoted by 
foil _ ( mt cosh y, k x , k y , m t sinh y). Let us assume that particles freeze out 
at a constant longitudinal proper-time ry , for the sake of simplicity. This im- 
plies freeze-out at a constant pressure, but at a space-time rapidity dependent 
temperature and density, and makes it possible to continue the calculation 
analytically. The source function of locally thermalized relativistically flowing 
particles in a Boltzmann approximation can be written as 

S(x,k) = C(rj) m t cosh(r] — y) n(x) exp (—k^u^/T) S(r — Tf), (31) 

where C(rj) is an 77 dependent normalization factor, given by the condition 
that / dk/E S(x, k) = n(x)5(r — r/), which implies that 

C(rj) = {47rm 2 T(r / ,7 ? )K 2 [m/T(r / ,r ? )]}^ 1 , (32) 

where K v (z) = f£°dzexp(— zcosht) cosh(z/t) is the modified Bessel function 
of the second kind. 

The single particle spectrum can be calculated from the emission function as 

d 3 N r 

E-— = / Tdrdr]S(x,k). (33) 



Substituting our family of new solutions, and using T(x) = 1/V(x), we obtain 
S(x,k) = C{rj) m t cosh(r] — y) n(x)/^(x, k) (34) 
/ B (a;,k)=exp r V( ^ ) 5{T-T f ). (35) 



to 
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We are interested in the coupling between the measurable rapidity distribution 
and the rapidity dependence of the effective temperature in the transverse 
directions as obtained from our new family of solutions. We assume that V(s) 
is a slowly varying function, i.e. d log V(s)/ds <C 1 in the region of interest. 
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This assumption implies that the point of maximal emissivity is located at 
fj = y with correction terms of 0(d log V(s) / ds) The measurable single-particle 
spectra can be written as 

E ^T = 2 C(y)n t V ( t -^y) K x [m t /T cS (y)], (36) 



dk yy> u u V Zl 

dN / tanh 



™ = 1IoioV . (37) 



where 



T C M= , \ a \ To I-) 1 "* ■ (38) 

y( ta gn V r // 

Note that the V function is a free fit function that describes the measurable 
rapidity distribution, including characteristic scales of the size of Z . 

We see that the slope parameter for transverse mass distribution T eff is related 
to the rapidity distribution as 

Figures 1 and 2 illustrate the calculated behavior of the effective temperature 
distribution as a function of rapidity for a single Gaussian-like and a double 
Gaussian-like ansatz for the measurable rapidity distribution. 

An interesting aspect of this new 1+1 dimensional solution is that the shapes 
of the rapidity distribution dN/dy and temperature distribution are coupled: 
the larger the rapidity density, the smaller the effective temperature. Choosing 
the effective temperature distribution T e s(y) to be flat, we recover the Hwa- 
Bjorken 1+1 dimensional solution, and the dN/dy rapidity distribution also 
becomes flat, rapidity independent. This behavior is expected to appear in 
high energy heavy ion collisions in the infinite bombarding energy limit. 



5 Cylindrically symmetric solutions 



In this section, we describe a new family of exact analytic solutions of rel- 
ativistic hydrodynamics, with cylindrically symmetric flow, overcoming the 
second of shortcoming of the well known boost-invariant Hwa-Bjorken solu- 
tion [9,10]. However, we do not address both shortcomings simultaneously 
yet. The physical motivation for this study is to consider the time evolution 
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dn/dy, T eff /MeV/ 
1000 




-3-2-10 1 2 3 

Fig. 1. Rapidity distribution dN/dy and effective temperature distribution T e ^{y) 
as a function of rapidity y, as obtained from a new family of solutions of 
(1+1) dimensional relativistic hydrodynamics. Here we use the scaling function 
V(s) = (1 — s)^ 1 / 4 ), using a scale parameter Zq = tanh(4), noto — 900 and 
7b(W r /) 1/K = 200 MeV , corresponding to a single maximum in the rapidity distri- 
bution dN/dy. The analytic expressions are given by eqs. (58, 60,37,38). 




-3-2-10 1 2 3 

Fig. 2. Same as Fig. 1 but utilizing a different form of the scaling function, 
V(s) = y/l + 1.6s 4 - 2.6s 8 , using a scale parameter Zq = 1, noio = 800 and 
To^o/rf) 1 ^ = 200 MeV, corresponding to a two-peaked rapidity distribution. 
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of central collisions in ultra-relativistic heavy-ion physics within the frame- 
work of an analytic approach. From now on, x M = (t,r x ,r y ,r z ) = (t,r) and 
fc" = (E, k x , ky, k z ) = (E, k) with E 2 - k 2 = m 2 . 

As we are primarily interested in the effects of finite transverse size and the de- 
velopment of transverse flow, we assume that the longitudinal flow component 
is boost-invariant, 

V*(*,0 = y. (40) 

We search for self-similar solutions, that are scale dependent in the tra nsverse 
directions, and depend only on the transverse radius variable r t = ^r 2 + r 2 
through the scaling variable 

_i_ ^,2 

• = ^, (41) 



and the longitudinal proper time r z = ^Jt 2 — r 2 and assume that, in the frame 
where v z — (longitudinal proper frame), the transverse motion corresponds 
to a Hubble type of self-similar transverse expansion, 

<^=M) r " * t -' t - ) = Wt" (42) 

where R = dR(r z )/dr z and hereafter we will designate by starred symbols 
the variables in the longitudinal proper frame. We assume that the scale R 
depends on time only through the longitudinal proper time, t z . 

In a relativistic notation, the above form may be parametrized as 

■u M = (cosh £ cosh £, sinh£— , sinh£— , sinh£cosh£), (43) 

r% r t 

tanh^ = ^p\r t = v* t = %v t , or cosh^ = - 1 = = 7* , (44) 
R ^> V 1 - R2s 

cosh^ = — = 7; • (45) 
t z 

The space-time rapidity rj is still defined by eq. (20). For a scaling longitudinal 
flow we obtain ( = rj . Using the above ansatz for the flow velocity distribution, 
we find that the continuity equation is solved by the form 

„(t,r x ,r tf ,r,)=no(^)^S( a ), (46) 

where G(s) is an arbitrary non- negative function of the scaling variable s and 
no , r z0 and i? are normalization constants. We use the convention n = 
n(t , 0, 0, 0), r z0 = T z (t , r z0 ) and R = R(t z0 ), where r z0 is such that, together 
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with to , satisfies eq. (40). This implies that Q(s = 0) = 1. The temperature 
equation, eq. (9), is solved by 

nU., V J = T G g^)"V W . (47) 

The constants of normalization are T = T(t , 0, 0, 0) and .F(O) = 1. We find 
that the solution is independent of the form of the function J-'(s), provided 
that F(s) > 0. 

Using a similar technique as in section 3, we obtain a transcendental equation 
for R 2 , and s. This equation has a particular solution if 

R = R Q = Const. (48) 

In this case, the acceleration of the radius parameter vanishes, R — 0, and the 
solution is R = R + Rq(t z — r z0 ). The relativistic Euler equation reduces to 

i 1 + + 3 ^ 2 ) = 2(1 - sR 2 ) [\o g g( S )T(s)}' , (49) 

where the lhs depends only on t z while the rhs is only a function of the 
variable s, hence both sides are constant. This implies that y — Rq ,thus 

Ro = Rqt z o ■ Thus the origin of the time axis (fixed by the assumption of 
the scaling longitudinal flow profile) coincides with the vanishing value of the 
transverse radius parameters. 

The solutions can be casted in a relatively simple form by introducing the 
proper time r, 

r = H ~ r " = \A 2 -rl~r$- r l ( 50 ) 



Using this natural variable we find that 
r 

v = -, or — — . 

f ' T 



(51) 



Thus the velocity field of our solution corresponds to the flow field of the 
spherically symmetric scaling solution and to the Hubble flow of the Universe. 
However, in the scaling solution the temperature and the pressure distributions 
are dependent only on the proper time r, while in our case both the density and 
the temperature distributions are generally dependent on the scale variable s 
in the transverse direction. 

As the solution is relativistic, and it is defined in the positive light-cone, given 
by t > 0, we obtain a constraint for the transverse coordinate, r t < r z . This 
together with the solution for the scale R, implies that the scaling variable 
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has to satisfy the constraint sR^ < 1, which corresponds to the limitation that 
the velocity of the fluid can not exceed the speed of light. 

By substituting R = Rqt z into the Euler equation, eq.(49), one obtains 

log [(1 - sR 2 )^ +1 ^Q(s)T(s)} = , (52) 

which gives, together with the condition ^(O)JF(O) = 1, 

G(s)T(s) = (1-R 2 s)- 2 ^ 1+1 / K \ (53) 

In this family of solutions, the scaling functions for the temperature and the 
density distribution are thus not independent. However, a constraint is given 
for their product, hence one of them can be chosen as an arbitrary positive 
function. For clarity, let us introduce new forms of the scaling functions as 

T(s)=Hs)H-Rls) 2/K , (54) 
V(s)=g(s)(l-Rls) 2 . (55) 

Then the constraint can be casted to the simple form of V(s)T(s) = 1. This 
construction for the scaling functions of the transverse density and temper- 
ature profiles coincides with the method, that we developed for the solution 
of the relativisitic hydrodynamical equations in the (1+1) dimensional prob- 
lem, but here the transverse flow has a two-dimensional distribution, so the 
exponents and the scaling variables had to be re-defined accordingly. 

Let us summarize our new family of solutions of the 1+3 dimensional rel- 
ativistic hydrodynamics for cylindrically symmetric systems by substituting 
the results to the density, temperature and pressure profiles. 



We obtain 



s 


A 






n(t,r) 


fr z0 
= n °{- 


p{t,r) 




T(t,r) 


= 7i (t 



-, for |r| < t, (56) 

t 

for r t <r z , (57) 



3 



V(s), (5* 

3+3/k 



3/k l 



v( s y 



(59) 
(60) 



where p = n T . Note that the scaling variable s is invariant for boosts in 
the longitudinal direction, and it is rotation-invariant in the transverse direc- 
tion, but s is not boost-invariant in the transverse directions. Hence we have 



13 



generated cylindrically symmetric, longitudinally boost invariant solutions of 
relativistic hydrodynamics. In the longitudinal direction, these solutions are 
homogeneous, boost-invariant and also scale- invariant. Due to this reason, the 
observable rapidity distribution is 

dN 

— — = const, (61) 
dy 

a flat distribution, corresponding to the ultra-relativistic nature of the solution 
in the longitudinal direction (where y = 0.5 \og[(E + k z )/(E — k z )\ is the 
rapidity of a particle with four-momentum (E,k) and dn/dy is the rapidity 
distribution of particle density). 

A new hydrodynamical solution is assigned to each non-negative function V(s), 
similarly to the cases of the non-relativistic solutions of ref. [18] and the 1+1 di- 
mensional relativistic solution of the previous section. Note that the solutions 
are valid also for massive particles, the form of the solution is independent of 
the value of the mass m. The form of solutions depends parameterically on k, 
that characterizes the equation of state. 

We have obtained new solutions of the (1+3) dimensional relativistic hydro- 
dynamical equations which describe a self-similar, streaming flow. In the case 
of R = 1 and V(s) = 1 we recover the spherically symmetric scale-invariant 
solution. This means that, in this limiting case, the pressure, the density and 
the temperature profiles depend only on the proper time r. In general case, 
however, our solution depends not only on the characteristic scale R but also 
on the arbitrary scaling function V(s). 



6 Summary 



We have found a new family of both 1+1 dimensional, longitudinally expand- 
ing, and 1+3 dimensional, cylindrically symmetric, adiabatic solutions of rel- 
ativistic hydrodynamics with conserved particle number. These families of 
solutions solves the continuity equation and the conservation of the energy - 
momentum tensor of a perfect fluid, assuming simple equations of state, given 
by Eqs.(4) and (5). The mass of the particles m and k = de/dp = l/c 2 s are 
free parameters of the solution. The well-known scale-invariant solution, has 
been obtained in the m = approximation. Interestingly, our generalizations 
resulted in additional freedom in the solution. 

In the new 1+1 dimensional hydro solutions, the flow field coincides with that 
of the Hwa-Bjorken solution. In principle, the shape of the measurable rapid- 
ity distribution, dN/dy plays the role of an arbitrary scaling function in our 
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solution, and we obtain that the effective temperature of the transverse mo- 
mentum distribution becomes rapidity dependent. Assuming that dN/dy is a 
slowly varying function of the rapidity y, we find that the effective temperature 
is proportional to the inverse of the rapidity distribution, T e d(y) oc (dN/dy)^ 1 . 

In 1+3 dimensions, even the flow velocity field deviates from Hwa-Bjorken so- 
lution. We find that the only exact solution in the considered class corresponds 
to a scaling 3-dimensional flow, similar to the Hubble flow of the Universe. Al- 
though the pressure distribution is only proper-time dependent, this pressure 
is a product of the local number density and the local temperature, hence one 
of these can be chosen in an arbitrary manner. 

The essential result of our paper is that we found a rich family of exact analytic 
solutions of relativistic hydrodynamics that contain both a longitudinal Hwa- 
Bjorken flow (that is frequently utilized in estimations of observables in high 
energy heavy ion collisions) and a relativistic transverse flow (whose existence 
is evident from the analysis of the single particle spectra at RHIC and SPS 
energies [23,24,25,26]). 

Acknowledgments: T. Cs. would like to thank L.P. Csernai, B. Lukacs and J. 
Zimanyi for inspiring discussions during the initial phase of this work, and 
to Y. Hama, G. Krein and S. S. Padula for their kind hospitality during his 
stay at USP and IFT, Sao Paulo, Brazil. This work has been supported by a 
NATO Science Fellowship (T. Cs.), by the OTKA grants T026435, T029158, 
T034269 and T038406 of Hungary, the NWO - OTKA grant N 25487 of The 
Netherlands and Hungary, and the grants FAPESP 00/04422-7, 99/09113-3, 
02/11344-8, PRONEX 41.96.0886.00, FAPERJ E-26/150.942/99, and CNPq, 
Brazil. 



References 



[1] L.P. Csernai, Introduction to Relativistic Heavy Ion Collisions, John Wiley and 
Sons, 1994. 



[2] D. Teaney, J. Lauret and E.V. Shuryak, nucl-th/01 10037 



vl. 



[3] T. Hirano, K. Morita, S. Muroya and C. Nonaka, nucl-th/0110009 vl. 



[4] P.F. Kolb, U. Heinz, P.Huovinen, K.J. Eskola and K. Tuminen, |hep^ph /0103234 
v3, Nucl. Phys. A696 (2001) 197-215. 

[5] C.E. Aguiar, T. Kodama, T. Osada, Y. Hama, J. Phys. G27 (2001) 75-94. 

[6] C.E. Aguiar, Y. Hama, T. Kodama, T. Osada, Nucl. Phys. A698 (2002) 639- 
642. 

[7] I.M. Khalatnikov, Zhur. Eksp. Teor. Fiz. 27 (1954) 529. 



15 



[8] L.D. Landau, Izv. Akad. Nauk SSSR 17 (1953) 51; in "Collected papers of L.D. 
Landau" (ed. D. Ter-Haar, Pergamon, Oxford, 1965) p. 569-585. 

[9] R.C. Hwa, Phys. Rev. D10 (1974) 2260; 

C.B. Chiu and K.-H. Wang, Phys. Rev. D12 (1975) 272; 
C.B. Chiu, E.C.G. Sudarshan and K.-H. Wang, Phys. Rev. D12 (1975) 902; 
M. I. Gorenstein, V. I. Zhdanov and Yu. M. Sinyukov, Sov. Phys. JETP 47 
(78) 435. 

K. Kajantie and L. D. McLerran, Phys. Lett. B 119 (1982) 203. 
K. Kajantie and L. D. McLerran, Nucl. Phys. B 214 (1983) 261. 



[10! 

[11 
[12 

[13! 
[14 

[15! 
[16 

[17! 

[18 
[19 

[20! 

[21 
[22 

[23! 

[24 

[25 

[26! 



J.D. Bjorken, Phys. Rev. D27 (1983) 140. 
T.S. Biro, Phys.Lett. B474 (2000) 21-26. 
T.S. Biro, Phys.Lett. B487 (2000) 133-139. 

J. Bondorf, S. Garpman and J. Zimanyi, Nucl. Phys. A296 (1978) 320. 

J.N. De, S.I. A. Garpman, D. Sperber, J. P. Bondorf and J. Zimanyi, Nucl. Phys. 
A305 (1978) 226. 

T. Csorgo, B. Lorstad and J. Zimanyi; Phys. Lett. B338 (1994) 134; 
nucl-th/9408022. 



J. Helgesson, T. Csorgo, M. Asakawa and B. Lorstad, Phys. Rev. C56 (1997) 
2626. 



P. Csizmadia, T. Csorgo and B. Lukacs, nucl-th/9805006 , Phys. Lett. B 443 
(1998) 21. 



T. Csorgo, [nucTth/9809011 



S.V. Akkelin, T. Csorgo, B. Lukacs, Yu.M. Sinyukov and M. Weiner, Phys. 
Lett. B 505 (2001), 64. 

T. Csorgo, S.V. Akkelin, Y. Hama, B. Lukacs and Yu.M. Sinyukov, 
hep-ph/0108067 ( accepted for publication in Phys. Rev. C. 

T. Csorgo, |a7Xiv:hep-ph/0111139| 



T. Csorgo and B. Lorstad, Phys. Rev. C54 (1996) 1390; T. Csorgo and B. 
Lorstad, Nucl. Phys. A590 (1995) 465c. 

W. Broniowski and W. Florkowski, Phys. Rev. Lett. 87 (2001) 272302 
[arXiv:nucl-th/0106050| . 

W. Florkowski and W. Broniowski. Acta Phys. Polon. B 33 (2002) 1629 
|arXiv:nucl-th /0203058] . 

W. Florkowski and W. Broniowski, arXiv:nucl-th/0212052l 

W. Broniowski, A. Baran and W. Florkowski, arXiv:nucl-th/0212053 



T. Csorgo and A. Ster, |arXiv:rmc l-th/0207016 



16 



